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Abstract 

We consider the two-variable interlace polynomial introduced by Arratia, 
Bollobas and Sorkin (2004). We develop graph transformations which allow 
us to derive point-to-point reductions for the interlace polynomial. Exploiting 
these reductions we obtain new results concerning the computational complex- 
ity of evaluating the interlace polynomial at a fixed point. Regarding exact 
evaluation, we prove that the interlace polynomial is #P-hard to evaluate at 
every point of the plane, except on one line, where it is trivially polynomial 
time computable, and four lines, where the complexity is still open. This solves 
a problem posed by Arratia, Bollobas and Sorkin (2004). In particular, three 
specializations of the two-variable interlace polynomial, the vertex-nullity in- 
terlace polynomial, the vertex-rank interlace polynomial and the independent 
set polynomial, are almost everywhere ^^P-hard to evaluate, too. For the inde- 
pendent set polynomial, our reductions allow us to prove that it is even hard 
to approximate at every point except at 0. 



1 Introduction 

The number of Euler circuits in specific graphs and their interlacings turned out 
to be a central issue in the solution of a problem related to DNA sequencing by 
hybridization |ABCSOO] . This led to the definition of a new graph polynomial, the 
one- variable interlace polynomial |ABS04a] . Further research on this polynomial 
inspired the definition of a two-variable interlace polynomial q{G]x,y) containing 
as special cases the following graph polynomials: qj^{G;y) = q{G;2,y) is the origi- 
nal one- variable interlace polynomial which was renamed to "vertex-nullity interlace 
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polynomial", qR{G; x) = q{G] x, 2) is the new "vertex-rank interlace polynomial" and 
I{G; x) = q{G; 1,1 + x) is the independent set polynomial [ ABS04bj . 

Although the interlace polynomial q{G]x,y) is a different object from the cele- 
brated Tutte polynomial (also known as dichromatic polynomial, see, for instance, 
|Tut84j ). they are also similar to each other. While the Tutte polynomial can be 
defined recursively by a deletion-contraction identity on edges, the interlace polyno- 
mial satisfies recurrence relations involving several operations on vertices (deletion, 
pivotization, complementation). 

Besides the deletion-contraction identity, the so called state expansion is a well- 
known way to define the Tutte polynomial. Here the similarity to the two-variable 
interlace polynomial is especially striking: while the interlace polynomial is defined 
as a sum over all vertex subsets of the graph using the rank of adjacency matrices 
(see (12.11) ). the state expansion of the Tutte polynomial can be interpreted as a sum 
over all edge subsets of the graph using the rank of incidence matrices (see (14. ip ) 
|ABS04bl Section 1]. 

References to further work on the interlace polynomial can be found in |ABS04b] 
and |EMSn7] . 

1.1 Previous work 

The aim of this paper is to explore the computational complexity of evaluating the 
two- variable interlace polynomial q{G;x,y). For the Tutte polynomial this problem 
was solved in |J"VW90j : Evaluating the Tutte polynomial is #P-hard at any alge- 
braical point of the plane, except on the hyperbola {x — l){y — 1) = 1 and at a 
few special points, where the Tutte polynomial can be evaluated in polynomial time. 
For the two- variable interlace polynomial q{G; x,y), only on a one-dimensional sub- 
set of the plane (on the lines x = 2 and x = 1) some results about the evaluation 
complexity are known. 

A connection between the vertex-nullity interlace polynomial and the Tutte poly- 
nomial of planar graphs |ABS04al End of Section 7], [EMS07t Theorem 3.1] shows 
that evaluating q is T^P-hard almost everywhere on the line x = 2 (Corollary 14. 4p . 

It has also been noticed that q{G; 1, 2) evaluates to the number of independent 
sets of G |ABS04b| Section 5], which is #P-hard to compute |Val79] . Recent work 
on the matching generating polynomial [AM07] implies that evaluating q is #P-hard 
almost everywhere on the line x = 1 (Corollary I4.10p . 

^The independent set polynomial of a graph G is defined as I{G;x) = X]j>o *(^' i)^"' ' where 
i{G; j) denotes the number of independent sets of cardinality j of G. 
^See Section [221 for a precise definition. 
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A key ingredient of |JVW90] is to apply graph transformations known as stretch- 
ing and thickening of edges. For the Tutte polynomial, these graph transformations 
allow us to reduce the evaluation at one point to the evaluation at another point. 
For the interlace polynomial no such graph transformations have been given so far. 

1.2 Our results 

We develop three graph transformations which are useful for the interlace polyno- 
mial: cloning of vertices and adding combs or cycles to the vertices. Applying these 
transformations allows us to reduce the evaluation of the interlace polynomial at 
some point to the evaluation of it at another point, see Theorem 13.31 Theorem 13.51 
and Theorem 13.71 We exploit this to obtain the following new results about the 
computational complexity of q{G;x,y). 

We prove that the two-variable interlace polynomial q{G;x,y) is T^P-hard to 
evaluate at almost every point of the plane. Theorem 14.111 see also Figure [H Even 
though there are some unknown (gray, in Figure [1]) lines left on the complexity 
map for q, this solves a challenge posed in [ABS04bl Section 5]. In particular we 
obtain the new result that evaluating the vertex-rank interlace polynomial qniG; x) 
is #P-hard at almost every point (Corollary I4.12p . Our techniques also give a new 
proof that the independent set polynomial is #P-hard to evaluate almost everywhere 
(Corollary OnD. 

Apart from these results on the computational complexity of evaluating the in- 
terlace polynomial exactly, we also show that the values of the independent set poly- 
nomial (which is the interlace polynomial q{G;x,y) on the line x = 1) are hard to 
approximate almost everywhere (Theorem 15.41) . 

2 Preliminaries 

2.1 Interlace Polynomials 

We consider undirected graphs without multiple edges but with self loops allowed. 
Let G = (V, E) be such a graph and A CV. By G[A] we denote {A, {e\e G -E, e C 
A}), the subgraph of G induced by A. The adjacency matrix of G is the symmetric 
n X n-matrix M = {rriij) over F2 = {0,1} with rriij = 1 iff {i,j} G E. The rank 
of this matrix is its rank over F2. Slightly abusing notation we write rk{G) for this 
rank. This allows us to define the two-variable interlace polynomial. 
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Definition 2.1 ( |ABS04b] ). Let G = {V,E) be an undirected graph. The interlace 
polynomial q{G] x, y) of G is defined as 

q{G;x,y)= ^ (a; - - (2.1) 

ACV 

In Section [3] we will introduce graph transformations which perform one and the 
same operation (cloning one single vertex, adding a comb or a cycle to one single 
vertex, resp.) on every vertex of a graph. Instead of relating the interlace polynomial 
of the original graph directly to the interlace polynomial of the transformed graph, 
we will analyze how, say, cloning one single vertex changes the interlace polynomial. 
To express this, we must be able to treat the vertex being cloned in a particular 
way, differently from the other vertices. This becomes possible using a multivariate 
version of the interlace polynomial, in which each vertex has its own variable. Once 
we can express the effect of cloning one vertex by an appropriate substitution of 
the vertex variable in the multivariate interlace polynomial, cloning all the vertices 
amounts to a simple substitution of all vertex variables and brings us back to a 
bivariate interlace polynomial. This procedure has been applied successfully to the 
Tutte polynomial |Sok05l[BM06] . 

We choose the following multivariate interlace polynomial, which is similar to the 
multivariate Tutte polynomial of Sokal [Sok05] and a specialization of the multivari- 
ate interlace polynomial defined by Courcelle |Cou07j . 

Definition 2.2. Let G = {V, E) he an undirected graph. For each v & V let x^ 
he an indeterminate. Writing xa for Y^veA^-"' define the following multivariate 
interlace polynomial: 

P{G;u,^) = 5^ 

ACV 

Suhstituting each x^ in P{G; u, x) hy x, we ohtain another hivariate interlace poly- 
nomial: 

ACV 

An easy calculation proves that q and P are closely related: 

Lemma 2.3. Let G he a graph. Then we have the polynomial identities q{G; x, y) = 
P{G;^,y-l) and P{G;u,x) = q{G;ux + l,x + l). □ 
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2.2 Evaluating Graph Polynomials 

Given ^, G Q we want to analyze the following computational problem: 
Input Graph G 
Output q{G]C,,v) 

This is what we mean by "evaluating the interlace polynomial q at the point 
As an abbreviation for this computational problem we write 

which should not be confused with the expression q{G; C,, v) denoting just a value in 
Q. Evaluating other graph polynomials such as P, q^, q^ and I is defined accordingly. 

If Pi and P2 are computational problems we use Pi P2 (Pi P2) to denote a 
polynomial time Turing reduction (polynomial time many-one reduction, resp.) from 
Pi to P2. For instance, Lemma [2.31 gives 

Corollary 2.4. For e Q, v 1, we have q{i, v) ^"^ P( JEj, v-1). For fi,^ eQ 
we have P{fi,^) q{^^ + 1,^+1). □ 

Here Q denotes some finite dimensional field extension Q C Q C M, which has 
a discrete representation. As \/2 will play an important role but we are not able to 
use arbitrary real numbers as the input for a Turing machine, we use Q instead of 
Q or M. We fix some Q for the rest of this paper. This construction is done in the 
spirit of Jaeger, Vertigan, and Welsh [JVWQOj who also propose to adjoin a finite 
number of points to Q in order to talk about the complexity at irrational points. To 
some extent, this is an ad hoc construction, but it is sufficient for this work. 

3 Graph Transformations for the Interlace Poly- 
nomial 

Now we describe our graph transformations, vertex cloning and adding combs or 
cycles to the vertices. The main results of this section are Theorem 13.31 Theorem 13.51 
and Theorem 13.71 which describe the effect of these graph transformations on the 
interlace polynomial. 

3.1 Cloning 

Cloning vertices in the graph yields our first graph transformation. 
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Cloning one vertex 



Let G = {V, E) be a graph. Let a E V he some vertex (the one which will be cloned) 
and N the set of neighbors of a, V = V \ {a} and M = V \ N. The graph G 
with a cloned, Gaa, is obtained out of G in the following way: Insert a new isolated 
vertex a'. Connect a' to all vertices in A^. If a does not have a self loop, we are done. 
Otherwise connect a and a' and insert a self loop at a'. Thus, adjacency matrices of 
the original (cloned, resp.) graph are 



B 





a 


N 


M 


a 


b 


1 





N 


1 


An 


A12 


M 





A21 


A22 



and 





a' 


a 


N 


M 


a' 


h 


b 


1 





a 


b 


b 


1 





N 


1 


1 


An 


A12 


M 








A21 


A22 



resp, (3.1) 



where 6 = 1 if a has a self loop and 6 = otherwise. As the first column of Baa equals 
its second column, as well as the first row equals the second row, we can remove the 
first row and the first column of Baa without changing the rank. This also holds 
when we consider the adjacency matrices of G[A] {Gaa[A], resp.) instead of G {Gaa 
resp.) for A <OV'. Thus we have for any ACV' 

rk{Gaa[A])=rk{G[A]), (3.2) 

rkiGaa[A U {a, a'}]) = rA;(G,JA U {a}]) = rkiGaa[A U {a'}]) = rkiG[A U {a}]). 

(3.3) 

Let X = ixv)y^v(Gaa) be a labeling of the vertices of Gaa by indeterminates. 
Define X to denote the following labeling of the vertices of G: := Xy for all 

V G V, Xa := (1 + Xa){l + Xa') — 1 = Xa + Xa' + XaXa'- Then wc havc 

Lemma 3.1. P{Gaa]u,yi) = P{G]u,ys^). 
Proof. On the one hand we have 

P(Gaa;M,x) 

= ^ Xa(m''''^*^'''"^^'^ a;^^i'-fc(Gaa[Au{a}]) 2;^,yrfc(Gaa[Au{a'}]) + a;^a;^,^''fc(Gaa[Au{a,a'}])^ 



ACV 



X^(m"^(^[^]) + {Xa + Xa' + XaXa')u 



rk{G[Avj{a}])\ 



by (S. 



ACV 



6 



On the other hand we have 

ACV 

□ 

Cloning all vertices 

Fix some k. Given a graph G, the graph Gk is obtained by cloning each vertex of 
G exactly k — 1 times. Note that the result of the cloning is independent of the 
order in which the different vertices are cloned. For a G V{G) let ai, . . . , a^. be the 
corresponding vertices in Gk- For a vertex labeling x of G^ we define the vertex 
labeling X of G by = (1 + XaJ(l + Xaj) ■ ■ ■ (1 + Xa^) — 1 for a G V{G). Applying 
Lemma 13.11 repeatedly we obtain 

Lemma 3.2. P(Gfc; u, x) = P(G; m, X) . □ 

Substitution of x^ by x for all vertices v gives 

Theorem 3.3. Let G he a graph and Gk he ohtained out of G hy cloning each vertex 
of G exactly k — 1 times. Then 

P{Gk;u,x) = P{G;u,{l + x)'' -1). (3.4) 

□ 

As we will use it in the proof of Theorem 14. IH we note the following identity for 
q, which can be derived from Theorem 13.31 using Lemma 12.31 

k 1 

q{Gk; X, y) = q{G- (x - 1)^^ + 1, y'). (3.5) 

Theorem 13.31 also implies the following reduction for the interlace polynomial, 
which is the foundation for our results in Section |H 

Proposition 3.4. Let B2 = {0, —1, —2} and x he an indeterminate. For /i G Q, ^ G 
Q \ i?2 we have P{fi,x) P{fi,^). (For any fi E Q, we write P{^,x) to denote 
the following computational prohlem: given a graph G compute P{G; fi,x), which is 
a polynomial in x with coefficients in Q.) 
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Proof. Let /i and ^ be given such that they fulfill the precondition of the proposition. 
Given a graph G =: Gi with n vertices, we build 6*2, Gs, • . • , Gn+i, where Gi is 
obtained out of G by cloning each vertex i — 1 times. This is possible in time 
polynomial in n. By Theorem 13.31 a call to an oracle for P{n, C.) with input Gi gives 
us P{G] /i, (1 + — 1) for i = 1, . . . , n + 1. The restriction on ^ guarantees that for 
2 = 1, 2, 3, . . . the expression (1 + ^)* — 1 evaluates to pairwise different values. Thus, 
for P{G; fi, x), which is a polynomial in x of degree < n, we have obtained the values 
at n + 1 distinct points. Using Lagrange interpolation we determine the coefficients 
ofP(G;/i,x). □ 



3.2 Adding Combs 



The comb transformation sometimes helps, when cloning has not the desired effect. 
Let G = (V, E) be a graph and a & V some vertex. Then we define the fc-comb of 
G at a as Ga,k = U {ai, . . . , ak}, E U {{a, ai}, . . . , {a, ak}}), with ai, . . . , being 
new vertices. 

Using similar arguments as with vertex cloning, adding combs to vertices yields 
a point-to-point reduction for the interlace polynomial, too. 

Theorem 3.5. Let G he a graph and Gk he ohtained out ofG hy performing a k-comh 
operation at every vertex. Then 



P{Gk] u, x) = p{k, u, x)l^(^)lp(G; u, x/p{k, u, x)), 
where p{k, u, x) = (1 + xY{xu^ + 1) 



(3.6) 



xu 



Proof. The adjacency matrices of the original graph G (the graph Ga,k with a /c-comb 
at a, resp.) are 





a 


V 


a 


b 


c 


V 




A 



and 







02 . 




a 


V 


ai 








1 




02 








1 




ak 








1 




a 


1 


1 . 


. 1 


h 


c 


v 










All 



resp. 



(3.7) 



with empty entries being zero. Consider A C V{Ga,k)- Let M := A fl {a, ai. 
By (13.71) . the rank of Ga,k is related to the rank for G in the following way: 
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• lia^M, then rk{Ga,k[^]) = rk{G[A \ M]). 

• If a e M and M n {oi, . . . , ^ 0, then rk{Ga,k[A]) = rk{G[A \M]) + 2: Let 
w.Lo.g. Oi e M. Consider the adjacency matrix of and the following 
operations on it, which leave the rank unchanged. Using the first column we 
remove all Is in the a-row, except the 1 in the first column. Using the first row 
we remove all Is in the a-column, except the 1 in the first row. The resulting 
matrix B is a {k + \V\) x [k + \V\) matrix with Is at positions (a, fli) and 
(ai, a), the submatrix of An induced by A \ M in the lower right corner and 
zeros everywhere else. Thus rk{B) = rk{G[A\ M]) + 2. 

• If M = {a}, then rk{Ga,k[A]) = rk{G[A]). 

Letting r{A) := rk{G[A]) and ra{A) := rk{G[A U {a}]) for A C V, we see that 
P{Ga,k]u,y.) equals 

^Xa(u'^^\ Xs + XaU^- Xs) 

ACV 9CSC{ai,...,ak} 0C5C{ai,...,afc} 

" V ' 

=:p{k,u,x) 

+ XaU''"^'^^^ 

Note that p{k, u, x) does only depend on Xa, , ■ ■ ■ , Xa^. , but not on x„ for any v e V . 
As we have 

ACV' 

we conclude 

P{Ga,k; u, x) = p{k, u, x)P(G'; u, X), 

where Xy = Xy for v & V and X^ — \ ■ 

We can perform a A;-comb operation at every a & V and call the result G^. 
Substituting x iov Xy,v & Gk, concludes the proof. □ 

This yields 

Proposition 3.6. Let p{k,u, x) = (1 + x)'' {xu"^ + 1) — xu"^ . Let k he a positive integer 
andii,^eQ. Ifp{k,ii,i)^Q, we have P{ii,^/p{k,ii,0) Pil^.O- □ 
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3.3 Adding Cycles 



Let G = (V, G) be a graph and a E V some vertex. Consider the graph Ga,k = 
{V U {l,2,...,k - 1},E U {{a, 1}, {a, k - 1}} U {{i - \ 1 < i < k}), with 
1, 2, . . . , — 1 being new vertices. We say that Ga,k has been obtained out of G by 
adding a k-cjde to a. 

Theorem 3.7. Let G be a graph and Gk be obtained out of G by adding a k-cycle to 
every vertex. Then P{Gk',u,x) = pk{u,x)P{G;u,qk{u,x)/pk{u,x)) fork = 3,4 with 
Pziu, x) = l + 2x + ?>x'^v? , qsi^u, x) = x + x^v? , p^i^u, x) = 1 + 3x + x"^ + 2x'^v? + x^u^ 
and q4{u,x) = + 2x^u'^ + x'^u'^. 

Proposition 3.8. P(0, 1) P(0, -1) and P{i2, -4) P(/i, -2) for every 12 e Q. 

Proof. The first reduction follows from Theorem 13.71 adding 3-cycles, the second 
adding 4-cycles. □ 

Proof of Theorem \3. 7[ We use the same idea as in the proof of Theorem l3.5[ Consider 
the case of a 3-cycle added at vertex a. Let V^' = \ {a}. The adjacency matrix of 
Ga,3 is 





1 


2 


a 


v 


1 




1 


1 




2 


1 




1 




a 


1 


1 


h 


c 


v 






T 


All 



with empty entries being zero. Adding the second row to the first row and the second 
column to the first column and subsequently the first row to the third row and the 
first column to the third column does not change the rank and gives 





1 


2 


a 


V 


1 




1 






2 


1 








a 






b 


c 


V 








An 



This shows that rk{GaAA] = rA;(G[A\ {1, 2}]) +2 for all A, {1,2, a} CAC V{Ga,3)- 
Using arguments similar to this one and the ones in the proof of Theorem 13.51 we 
find that 

• rk{Ga,3[A]) = rk{G[A nV'] + 2 for all A, {a} C A C V{Ga,3) and either 1 e A 
or 2 e A, 
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rk{GaA^]) = rk{G[A] for all A, {a} C A C V{Ga,3) and {1, 2} n A = 0, 



. rk{Ga,3[A]) = rk{G[A n V']) + rkiP^iA n ^(Ps)]) for all A C V{Ga,3), a ^ A, 
where P2 is the the path with two vertices 1,2. 

Letting again r{A) := rk{G[A]) and r„(yl) := rk{G[AL} {a}]) for A C 1^', we see that 
P{Ga,3] u, x) equals 

= :(j3{«,x) 

which equals P3(m, x)P(G'; m, X) if we define X by = for t> G and = 
g3(M, x)/p3(m, x). We can use this identity for every vertex a and substitute Xa, 
a G V^, by a single variable x. This gives the statement of the theorem concerning 
3-cycles. For 4-cycles we proceed in a similar fashion. □ 



4 Complexity of evaluating the Interlace Polyno- 
mial exactly 

The goal of this section is to uncover the complexity maps for P and q as indicated 
in Figure [H While the left hand side (complexity map for P) is intended to follow 
the arguments which prove the hardness, the right hand side (complexity map for q) 
focuses on presenting the results. 

Remark 4.1. P(/i, 0) and P(l,^) are trivially solvable in polynomial time for any 
/X, e e Q, as P{G; /i, 0) = 1 and P{G; 1, = (1 + O'^' ■ □ 

Thus, on the thick black lines x = and m = 1 in the left half of Figure [Tj P can 
be evaluated in polynomial time. By Lemma [2.31 these lines in the complexity map 
for P correspond to the point (1, 1) and the line x = y, resp., in the complexity map 
for q, see the right half of Figure [H 



4.1 Identifying hard points 

We want to establish Corollary 14.41 and Remark |4. 51 which tell us, that P is ^P-hard 
to evaluate almost everywhere on the dashed hyperbola in Figure [1] and at (0, 1). To 
this end we collect known hardness results about the interlace polynomial. 
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in P 



Figure 1: Complexity of the interlace polynomials P and q. a = a/2, (3 = l/\/2 



Let t{G] X, y) denote the Tutte polynomial of an undirected graph G = {V, E). It 
may be defined by its state expansion as 

t(G;x,y)= (x-ir(^)"^(^)(y-l)l^l-^'(^\ (4.1) 

B<ZE(G) 

where r{B) is the F2-rank of the incidence matrix of G[B] = {V, B), the subgraph of 
G induced by B. (Note that r{B) equals the number of vertices of G[B] minus the 
number of components of G[-B], which is the rank of B in the cycle matriod of G.) 
For details about the Tutte polynomial we refer to standard literature |Tut84[ IBU921 
IWel93] ■ The complexity of the Tutte polynomial has been studied extensively. In 
particular, the following result is known. 

Theorem 4.2 ( [Ver05] ) . Evaluating the Tutte polynomial of planar graphs at (^,^) 
IS i^P-hard for all ^ e Q except for ^ G {0, 1, 2, 1 ± ^2}. 

We will profit from this by a connection between the interlace polynomial and 
the Tutte polynomial of planar graphs. This connection is established via medial 
graphs. For any planar graph G one can build the oriented medial graph Gm, find 
an Euler circuit C in Gm and obtain the circle graph H of G. The whole procedure 
can be performed in polynomial time. For details we refer to |EMS07j . We will use 
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Theorem 4.3 ( |ABSn4a[ End of Section 7]; |EMSn7[ Theorem 3.1]). Let G be a 

planar graph, Gm he the oriented medial graph of G and H he the circle graph of 
some Euler circuit C of Gm- Then q{H; 2, y) = t{G] y, y). Thus we have t{v, v) 
P(^j,v — 1), where t{v,v) denotes the prohlem of evaluating the Tutte polynomial 
of a planar graph at {v,v). 

Proof. See the references for q{H; 2, y) = t{G; y, y) and use Lemma [231 □ 

We set a = 72 and (3 = 1/V2. Let Bi = {±l,±/3,0}. Theorem and 
Theorem 14.31 yield 

Corollary 4.4. Evaluating the vertex-nullity interlace polynomial qn is ^P-hard 
almost everywhere. In particular, we have: 

• The prohlem QnC^) is trivially solvable in polynomial time. 

• For any v E Q \ {0, 1,2, 1 ± a} the problem qN{v) = q{'2,v) is ^P-hard. Or, 
in other words, for any /U G Q \ -Bi the problem P{^^ V/^) i^P-hard. □ 

Remark 4.5. P(0, 1) is ^P-hard, as P(G; 0, 1) equals the number of independent 
sets ofG, which is jj^P-hard to compute lVal79^ . □ 



4.2 Reducing to hard points 

The cloning reduction allows us to spread the collected hardness over almost the 
whole plane: Combining Corollary 14.41 and Remark 14.51 with Proposition 13.41 we 
obtain 

Proposition 4.6. Let Bi = {±l,±/9,0} and B2 = {0,-1,-2} (as defined on 
Pages {IE and\^ resp.). Let e ((Q \ Bi) U {0}) x (Q \ B2). Then Pifi,0 

is ^P-hard. □ 

This tells us that P is T^P-hard to evaluate at every point in the left half of 
Figure [T] not lying on one of the seven thick lines (three of which are solid gray ones, 
two of which are solid black ones, and two of which are dashed brown ones). Using 
the comb and cycle reductions we are able to reveal the hardness of the interlace 
polynomial P on the lines x = ~1 and x = —2: 

Proposition 4.7. For G (Q \ -Bi) U {0} the problem P{lJ., —1) is ^P-hard. 



13 



Proof. For yU = we use Proposition 13.81 and Remark 14.51 If /i 7^ 0, we can use 
Proposition 13.61 which yields P(yU, — P(/i, — 1). For /i = ±1 this reduces 
(±1,-1) to itself. For fi = ±P this reduces -2) to -1) and -2) to 
(— /3, —1). For other fi this gives a reduction of some point, which is already known 
as #P-hard by Proposition 14. 6[ to (/i, —1). □ 

Proposition 4.8. For jj, e {Q\ Bi) U {0} the problem P(/i, —2) is ^P-hard. 

Proof. Use Proposition 13.81 and Proposition 14.61 □ 

4.3 Summing up 

First we summarize our knowledge about P. 
Theorem 4.9. Let [3 = I/V2. 

1. P(/i, ^) is computable in polynomial time on the lines yU = 1 and ^ = 0. 

2. For (/i, e (Q \ {-1, /5, 1}) X (Q \ {0}) the problem P{fx, is i^P-hard. 

Proof. Summary of Remark 14.11 Proposition 14.61 Proposition 14.71 Proposition 14.81 
The hardness of P(0, — 1) follows from Corollary 14.101 □ 

In particular we obtain the following corollary about the complexity of the inde- 
pendent set polynomial, which also follows from jAM07j . 

Corollary 4.10. Evaluating the independent set polynomial I {X) = P(0, A) = g(l, 1+ 
A) is ^P-hard at all X E Q except at X = 0, where it is computable in polynomial 
time. 

Now we turn to the complexity of q, see also the right half of Figure [B 

Theorem 4.11. The two-variable interlace polynomial q is ^P-hard to evaluate 
almost everywhere. In particular, we have: 

1. q{C,,v) is computable in polynomial time on the line = v. 

2. Let ^ G Q\{1} and x be an indeterminate. Then 1) is as hard as computing 
the whole polynomial q{x, 1). 

3. q{C,,v) is ^P -hard for all 

(e, v) e {{^,v) eQ"^ \vj^ + 1 and ±V2(^ - 1) + 1 and v 1}. 
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Proof of Theorem 4-1 A (Sketch). ([T]) and ([3]) follow from Remark |4. H and Theorem l4.9l 



using Lemma [231 For ^ 7^ 1, (13. 5p gives q{Gk]^, 1) = q{G; k{C, — 1) + 1, 1), which 
yields enough points for interpolation in the same way as in Proposition 13.41 using 
/c = 1,2,3, ... This proves p. □ 

Theorem 14.111 implies 

Corollary 4.12. Let (3 = l/\/2. Evaluating the vertex-rank interlace polynomial 
qniG; x) is jj^P-hard at all & Q except at ^ = 0,1 — P,l + P (complexity open) and 
^ = 2 (computable in polynomial time). □ 



5 Inapproximability of the Independent Set Poly- 
nomial 

Provided we can evaluate the independent set polynomial at some fixed point, vertex 
cloning (adding combs, resp.) allows us to evaluate it at very large points. In this 
section we exploit this to prove that the independent set polynomial is hard to 
approximate. Similar results are shown in jGJ07j for the Tutte polynomial. 

Definition 5.1. Let A G Q and e > 0. By a randomized 2"^ " -approximation algo- 
rithm for /(A) we mean a randomized algorithm, that, given a graph G with n nodes, 
runs in time polynomial in n and returns I{G; A) G Q such that 

Pr[2-"'"'J(G; A) < J(G; A) < 2"'"'J(G; A)] > ^. 

In |GJ07j . (non)approximability in the weaker sense of (not) admitting an FPRAS 
is considered. 

Definition 5.2. Let A G Q. ^4 fully polynomial randomized approximation scheme 
(FPRAS) for /(A) is a randomized algorithm, that given a graph G with n nodes and 
an error tolerance e,0 < e < 1, runs in time polynomial in n and 1/e and returns 
i{G; A) G Q such that 

Pr[2-"/(G; A) < /(G; A) < 2'I{G; A)] > ^. 

Lemma 5.3. For every A G Q, 7^ |1 + A| 7^ 1, and every e, < e < 1, there is no 
randomized polynomial time 2"' ^-approximation algorithm for I {X) unless RP = NR. 
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Theorem 5.4. For every A G Q\{0} and every e, < e < 1, there is no randomized 
polynomial time 2" " -approximation algorithm (and thus also no FPRAS) for /(A) 
unless RP = NP. 



Proof. Lemma [5.31 gives the inapproximabihty at A G Q \ {—2, —1,0}. By (13. 6p we 
could turn an approximation algorithm for /(— 2) into an approximation algorithm for 
/(2) which would imply RP = NP by Lemma 15731 For /(—I) we use Theorem 13. 71 □ 



Proof of Lemma [JTM Fix A G Q, 7^ |1 + A| 7^ 1, and £, < e < 1. Assume we 
have a randomized 2" ^-approximation algorithm A for /(A). Given a graph G, 
Theorem 13.31 and Theorem 13.51 resp., will allow us to evaluate the independent set 
polynomial at a point ^ with |^| that large, that an approximation of /(G; C) can be 
used to recover the degree of I{G;x), which is the size of a maximum independent 
set of G. As computing this number is NP-hard, a randomized 2" "-approximation 
algorithm for I{G; A) would yield an RP-algorithm for an NP-hard problem, which 
implies RP = NP. 

Let G = {V, E) be a graph with \ V\ = n. We distinguish two cases. If 1 1 + A| > 1, 
we choose a positive integer / such that (n/)^^^ > and with ^ := (1 + A)' — 1 we 
have 

As A and e are constant, this can be achieved by choosing / = poly(n). If < 
|1 + A| < 1, we choose a positive integer / such that with ^ := -jj^^ (15. ip holds. 

By Theorem (Theorem [331 resp.) we have /(G;0 = HGi;\) (/(G; = (1 + 
A)^''^'/(Gi; A), resp.). Algorithm A returns on input Gi within time poly(n/) = 
poly(n) an approximation /(G;;A), such that with I{G;^) := I{Gi;\) (/(G;^) := 

' ^^^P-) ^^^^ 



7(G;0</(G;0<2('^')"^/(G;0 (5.2) 

with high probability. 

Let c be the size of a maximum independent set of G, and let N be the number 
of independent sets of maximum size. We have 



/(G; x) = Nx" + 



o<i<c-i 



i{G-3)x^ 



and thus 



-N 



0<j<c-l 



< c2"|^|-^ < 2 



jlogn.+?i|^|-l 



(EU 1 

< -. 



(5.3) 
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If we could evaluate I{G] ^) exactly, we could try all c G {1, . . . , ra} to find the one for 
which ^^^^^ is a good estimation for A^, 1 < < 2". This c is unique as |,^| > 2"^ 
The following calculation shows that this is also possible using the approximation 
algorithm A. 

Using A we compute N{c) := ^^^p^ for all c G {1, ... , n}. We claim that c is the 
unique c with 

2-(n0i---i < < 2{«0^-^+n+i. (5.4) 
Let us prove this claim. As 1 < < 2" and by (15. 3p . we know that 

I < < 2"+^ (5.5) 

Thus, by (El, c = c fulfills 

On the other hand, when c < c — 1 we have 

\N{d)\ > 2-("') > 2("') +"+\ 

When c > c + 1 we have |iV(c)| < 2^'-"'^^ by similar arguments. This shows 
that any integer c,Cy^c, does not fulfill (15.41) . Thus, c can be found in randomized 
polynomial time using A. □ 
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